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ABSTRACT 

Context. Since 1998, a planet-search program around main sequence stars within 50 pc in the southern hemisphere has been carried 
out with the CORALIE echelle spectrograph at La Silla Observatory. 

Aims. With an observing time span of more than 14 years, the CORALIE survey is now able to unveil Jovian planets on Jupiter's 
period domain. This growing period-interval coverage is important for building formation and migration models since observational 
constraints are still weak for periods beyond the ice line. 

/Wefhods. Long-term precise Doppler measurements with the CORALIE echelle spectrograph, together with a few additional ob- 
servations made with the HARPS spectrograph on the ESO 3.6 m telescope, reveal radial velocity signatures of massive planetary 
companions on long-period orbits. 

Results. In this paper we present seven new planets orbiting HD 27631, HD 98649, HD 106515A, HD 166724, HD 196067, 
HD 219077, and HD 220689, together with the CORALIE orbital parameters for three already known planets around HD 10647, 
HD 30562, and HD 86226. The period range of the new planetary companions goes from 2200 to 5500 days and covers a mass do- 
main between 1 and 10.5 Mj^p. Surprisingly, five of them present very high eccentricities above e > 0.57. A pumping scenario by 
Kozai mechanism may be invoked for HD 1065 15Ab and HD 196067b, which are both orbiting stars in multiple systems. Since the 
presence of a third massive body cannot be inferred from the data of HD 98649b, HD 166724b, and HD 219077b, the origin of the 
eccentricity of these systems remains unknown. Except for HD 10647b, no constraint on the upper mass of the planets is provided by 
Hipparcos astrometric data. Finally, the hosts of these long period planets show no metallicity excess. 

Key words, stars: planetary systems - stars: binaries: visual - techniques: radial velocities - Stars: individual: HD 27631 - Stars: 
individual: HD 98649 - Stars: individual: HD 1065 15A - Stars: individual: HD 166724 ~ Stars: individual: HD 196067 - Stars: indi- 
vidual: HD 219077 - Stars: individual: HD 220689 - Stars: individual: HD 10647 - Stars: individual: HD 30562 - Stars: individual: 
HD 86226 



1. Introduction 



Send offprint requests to: Maxime Marmier, e-mail: 
Maxime . MarmierOunige . ch 

* The CORALIE radial velocity measurements discussed in this 
paper are only available in electronic form at the CDS via anony- 
mous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u- 
strasbg.fr/cgi-bin/qcat ?J/A+ A/ 

** Based on observations collected with the CORALIE echelle spec- 
trograph mounted on the 1.2m Swiss telescope at La Silla Observatory 
and with the HARPS spectrograph on the ESO 3.6 m telescope at La 
Silla (ESO, Chile) 



The CORALIE planet-search program in the southern hemi- 
sphere ( Udry et al.,20 00) has been going on for more than 14 
years (since June 1998). This high-precision radial velocity sur- 
vey makes use of the CORALIE fiber-fed echelle spectrograph 
mounted on the 1.2 meter Euler Swiss telescope located at La 
Silla Observatory (ESO, Chile). The stars being monitored are 
part of a volume-limited sample. It contains 1647 main sequence 
stars from F8 down to KO within 50 pc and has a color-dependent 
distance limit for later type stars down to MO with the faintest 
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target not exceeding V=10. Today the survey comprises more 
than 35 000 radial- velocity measurements on the 1647 stars of 
the program. 

If additional constraints are added to reject the binaries, the high- 
rotation stars and the targets with a chromospheric activity index 
greater than log 7?^^ = -4.75 (or more than -4.70 for K stars), 
half of the original sample remain (822 stars), constituting the 
best candidates for a Doppler planet search survey. This sub- 
sample is observed with a precision of 6 ms" ' or better, while the 
remaining part of the sample suffers from a degraded accuracy 
due to stellar jitter and/or the widening of the cross-correlation 
function (induced by high stellar rotation velocities). However, 
ninety percent of the targets are still monitored with a long-term 
precision better than 10 ms"' . The statistical distribution of num- 
ber of measurements per star for the entire sample and the sub- 
sample (of the best candidates for Doppler-planet search survey) 
are represented in Fig.[T] Even if most of the discovered planets 
were found around stars belonging to the "quiet" subsample (for 
more details see Mayor et al.||201 1) , the remaining part of the 
survey has also been observed well. The median of the distribu- 
tion for the entire sample and the subsample are respectively 15 
and 17 Doppler measurements, while ninety percent of the stars 
have at least six data points. 

So far CORALIE has allowed (or has contributed to) detection 
of more than eighty extrasolar planet candidates. In this paper 
we report the discovery of seven long-period and massive plan- 
ets orbiting HD 27631, HD 98649, HD 1065 15 A, HD 166724, 
HD 196067, HD 219077, and HD 220689. We also provide up- 



dates for three known systems around HD 10647 (Butler et al. 
2006), HD 30562 ( [Fischer et al.l2009| , and HD 86226 (MmiiitT 
et al. 2009). After 14 years of activity, the CORALIE survey 



is now able to unveil planetary candidates with long orbital 
periods. HD 106515Ab and HD 196067b, with periods around 
ten years, have just completed one full revolution since the 
beginning of our observations. In the case of the companions 
orbiting HD 98649, HD 166724, and HD 219077, with periods 
around five thousand days, the orbital phase coverage is not 
complete yet. This growing period-interval coverage is very im- 
portant with regard to formation and migration models since 
only 28 of the 659 exoplanets (including 5 systems presented 
in this paper) with known radial velocity or/and astrometric or- 
bits have periods longer than nine years. Interestingly the ec- 
centricities of these five systems are all above 0.57. With eccen- 
tricites of 0.85, 0.734, and 0.770, HD 98649b, HD 166724b, and 
HD 219077b, respectively, are the three planetary companions 
with the most eccentric orbits on periods longer than five years 
(Fig.fTSll. Finally, with a mass higher than 6Mjup, HD 98649, 



HD 1(565 15 A, HD 196067, and HD 219077 belong to the upper 
15% of the planetary mass distribution, and they enrich the sta- 
tistical sample towards the region of the brown-dwarf desert. 
The paper is organized as follows. The new host stars properties 
are summarized in Sect.|2] Section[3]presents our radial-velocity 
data and the inferred orbital solution of the newly detected com- 
panions. Section|4]presents the stellar properties and CORALIE 
updated parameters of the three already known planets. These 
results are discussed in Sect.|5]with some concluding remarks. 

2. Stellar characteristics 



Spectral types and photometric parameters were taken from 
HIPPARCOS [ESA| ([1997]). The astrometric parallaxes come 
from the improved Hipparcos reduction (van Leeuwen||2007 1. 
The atmospheric parameters Teff, lo g g, [Fe/H] are from |Santos 
|et al.| ( |200T] l and |Sousa et al.| ( (2008l ) while the wsin (i) was com- 





# measurements 



# measurements per star 



Fig. 1. Statistical distribution of the number of measurements 
in the CORALIE volume-limited sample (1647 stars). The blue 
curve represents the subsample (822 stars) of the best candidates 
for a Doppler-planet search survey (low chromospheric activity 
index, no high-rotation stars, and no binaries), while the black 
line is for the entire sample. The histograms of the number of 
Doppler measurements are plotted on the right, and the curves 
in the left panel are the inverse cumulative distribution functions 
of measurements (i.e. the number of stars with at least N data 
points). The dashed lines represent the median of the distribu- 
tions for all them and for the subsample (respectively 15 and 17 
Doppler measurements). 



puted using the Santos et al. ( 2002 1 calibration of CORALIE's 



cross-correlation function (CCF). Using the stellar activity in- 
dex log R'ffj^, derived from HARPS or CORALIE high signal- 
to-noise spectra, we derived the rotational period from the em- 
pirical calibration described in ,Mamajek & Hillenbrand (2008} . 



The bolometric correction was computed from Flower ( 1996^ 



using the spectroscopic Tgff determinations. A Bayesian esti- 
mation method describ ed in|da Silva et al.| ( |2006| ) and theoreti- 
cal isochrones from Girardi et al.'(2000) were used to estimate 
stellar ages, masses and radii{^The observed and inferred stellar 
parameters are summarized in Table[T] 



2.1. HD 27631 (HIP 20199) 

HD 27631 is identified as a G3IV/V in the Hipparcos cata- 
log [ESA] ( [1997) with a parallax of tt = 22.0 + 0.7 mas and 
an apparent magnitude V - 8.26. The spectroscopic analysis 
of the CORALIE spectra results in an effective temperature of 
T,ff = 5737 + 36, a stellar metallicity of [Fe///] = -0.12±0.05, 
and a surface gravity \ogg - 4.48±0.09. Using these parameters 
and theoretical isochrones from [Girardi et al. (2000), we obtain 
a mass of M, = 0.94 + 0.04 Mq with an unconstrained age of 
4.4 + 3.6 Gyr No emission peak is visible for HD 27631 at the 
center of the Ca II K absorption line region (Fig.[2]). We derive a 
lo g^Hr-~^-94 ± 0.08 compatible with the value of -4.91 found 
by|Henry et"aL] ([1996)1. 



2.2. HD98649 (HIP55409, LTT4199) 

HD 98649 is a chromospherically-quiet G4 dwarf with an astro- 
metric parallax of tt = 24.1+0.8 mas, which sets the star at a dis- 
tance of 41 .5 ± 1 .4 pc from the Sun. The apparent magnitude V = 
8.00 implies an absolute magnitude of M^, - 4.91. According to 
the Hipparcos catalog, the color index for HD 98649 is B-V = 
0.658. Using a bolometric correction BC - -0.083 and the solar 



' The web interface (PARAM 1.1) for the Bayesian estimation of 
stellar parameters can be found at htt p://stev.oapd.inaf. it/cgi-bin/param| 
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Table 1. Observed and inferred stellar parameters for HD 27631, HD 98649, HD 106515A, HD 166724, HD 196067, HD219077, 
and HD 220689. 
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Parameter from HIPPARCOS ('ESA |1997| ). '*' Parallax from the nevy Hipparcos reduction |van Leeu wenl20d7|. 

~ . .. ^ HA RPS spectra ( |Sousa et al.|2008| . <"' Parameter derived using CO KALlb (J(Jh (Santosj 

aDf2000l models. '^' The bolometric correction is computed from Flower ( 1996). ''" Parameter 



Notes. "" 

using COR ALIE spectra 



2001 1. '^' Parameter deri ved from j Girardi et 



derived using HARPS spectra i Lovis et"arpOTT|~^From the calibration of the rotational period vs. activity ^Mamajek & Hillenbrand_2008^. 



absolute magnitude M^o/ - 4.83 (Lejeun e et al.|[T998] l, we ob- 
tain a luminosity L - 0.86 Lq. Our spectral analysis results in an 
effective temperature of Tgff - 5759 ± 35 and a stellar metal- 
licity of \FelH\ - -0.02 ± 0.03. We finally derived a mass of 
= 1.00 ± 0.03 Mq and an age of 2.3 + 2.0 Gyr HD 98649 
has an activity index l og /^'^^ - -4.9 9 ± 0.08, in agreement with 
the value published by [Jenkins et al.| ( |2008| l. The projected rota- 
tional velocity is vi/n (;) = 1.19 kms'' with a rotational period 
of Pro/ = 27 ±4 days. 



2.3. HD1 0651 5A (HIP 59743, CCDMJ1 21 51-71 6A, 
LTT4599) 

HD 1065 15A is the primary star of a triple visual system located 
in the Virgo constellation at 35.2 + 1.2 pc from the Sun. It 
has spectral type G5, an apparent magnitude V - 7.35 with 
an effective temperature of Te/f = 5362 + 29 and a stellar 
metallicity of [Fe/H] - 0.03 + 0.02. The derived mass is 
M, = 0.97 + 0.01 Mo with an age of 11. 7 + 0.2 Gyr The star is 
thus a very slow rotator \sin(i) < 1 kms"' with a low activity 
index log/?^^=-5.09+0.09. These values agree closely with the 
result published by |Desidera et al. (2006 1 (vsin (i) = 0.6 kms 



-5.04). HD106515A's first companion LTT4598 
(CCDMJ12151-716B) is a G8 dwarf 7.5 arcsecond away with 
an apparent magnitude V=8.3. 

Our spectroscopic analysis results in Tgff - 5324 ± 43, 
[Fe/H] = 0.09 ± 0.03, log^ = 4.47 ± 0.08, and a derived mass 
M» = 0.89 Mq. Owing to their small angular separation, the 



two stars have the same name in the Hipparcos catalog. Gould 



& Chaname (2004i derived an Hipparcos-based parallax for 
LTT4598 and confirm that the two stars have common proper 
motion, implying that the pair is gravitationally bound. Using 
the positions given by Gould & Chaname (2004), we obtain 
a projected binary separation of 310 AU, which leads to an 
estimated semi-major axis of 390 AU with the statistical relation 
a/r = 1.26 ( [Fischer & Marcy]|1992| ). A second companion at 
a separation of 98.7 arcsecond with a magnitude V = 10.3, 
identified as CCDM J12151-716C (BD- 06 3533, TYC4946- 
202-1) is Hsted in the CCDM catalog ( [Dommanget & Nys| 



'2002). However, the Tycho-2 catalog ('H0g et al.'2000) provides 
different proper motion than for HD 1065 15 A. The pair is thus 
probably not physical. 



2.4. HD 166724 (HIP 89354) 

HD 166724 has a spectral type KOIV/V in the Hipparcos cat- 
alog with an apparent "V band magnitude V - 9.33 and a 
parallax n - 23.6 + 1.2 mas. The spectral analysis results in 
an effective temperature of Tg/f = 5127 + 52 and a stellar 
metallicity of [Fe/H] - -0.09 ± 0.03. We derived a mass of 
M« = 0.81 ± 0.02 Mo with a luminosity of L = 0.31 Lo and an 
unconstrained age of 4.0 + 3.8 Gyr. HD 166724 is slightly ac- 
tive with a mean stellar activity index derived form the HARPS 
spectra of log R'^j^ - -4.73 + 0.03. A strong re-emission peak is 
visible for HD 166724 at the center of the Ca II K absorption line 
(Fig|2]i, which indicates the significant chromospherical activity 
of the star 



2.5. HD 196067 (HIP 1021 25, CCDM J2041 7-7521 A, 
LTT8159) 

HD 196067 is the primary component of a bright visual binary 
located at 44 + 5 pc from the Sun in the Octans constellation. Its 
spectral type is Gl with an apparent magnitude V = 6.51. The 
effective temperature and metallicity are T^ff - 6017 ± 46 and 
[Fe/H] = 0.18 + 0.04. We derived a mass M, = 1.29 + 0.08 Mo, 
a radius of P = 1.73 ± 0.21 R^, and an age of 3.3 + 0.6 Gyr. The 
star is chromospheiically inactive with a log - -5.11 +0.08, 
leading to a rotational period of 26 + 3 days. HD 196067's visual 
companion is HD 196068 (HIP 102128), a Gl dwarf with an ap- 
parent magnitude V=7.18. The apparent separation of the pair is 
16.8 arcsecond, which explains the two entries in the Hipparcos 
catalog. However, the two stars have very similar proper mo- 
tions and are known to be gravitationally bound (see e.g. |Gould| 
& Chaname 2004). This conclusion is strengthened by the fact 
that the CORALIE average radial- velocities of the two stars are 
identical within uncertainties. The Hipparcos positions lead to 
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[Fe/H] = -0.13 + 0.01, in agreement with the values published 
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Fig. 2. Call K absorption Hne region {A = 3933.6 A) of the 
co-added CORALIE spectra. This region is contaminated by 
a few thin thorium emission lines that were filtered while co- 
adding the spectra. No significant re-emission lines are present 
at the bottom of the Call absorption region for the stars dis- 
played in the upper panel, from top to bottom HD 27631, 
HD 98649, HD 106515A, HD 196067, HD 219077, HD 220689, 
HD 10647, HD 30562, and HD 86226 . The large re-emission 
for HD 166724 (lower panel) indicates the significant chromo- 
spheric activity of the star. 



a projected binary separation of 740 AU. This translates into 
an estimated binary semi-major axis of 932 AU. The spectro- 
scopic analysis of HD 196068 results in T,ff = 5997 + 44, 
[Fe/H] = 0.34 ± 0.04, logg = 4.38 ± 0.05, and vsin (/) = 1.21. 
These parameters lead to a stellar mass of - 1.18 + 0.07 M© 
with a poorly constrained age of 2.4 ±1.6 Gyr, and a radius of 
R = 1.34 ± 0.23 /?o. 



2.6. HD 21 9077 (HIP 11 4699) 

HD 219077 is a bright G8 dwarf with a magnitude V = 6.12 
and an astrometric parallax of tt = 34.07 ± 0.37, which cor- 
responds to an absolute magnitude of M,, = 3.78 and a lumi- 
nosity of 2.66 Lq. The star is chromospherically quiet with a 
mean stellar activity index of log/?^^ - -5.12 ± 0.01, lead- 
ing to a calibrated rotational period of 49 ± 4 days. The spec- 
tral analysis made by Sousa et al. ( 2008 1 results in an effective 
temperature of T^ff = 5362 + 18 and a stellar metallicity of 



by Ghezzi et al. (2010 1. Using the theoretical isochrones, we de- 
rived a mass of M. = 1.05 ± 0.02 Mq and an age of 8.9 + 0.3 
Gyr. 



2. 7. HD 220689 (HIP 1 15662) 

HD 220689 has a spectral type G3 V with an astrometric parallax 
of ;r = 22.4 ± 0.7 mas, which sets the star at a distance of 44.6 ± 
1 .4 pc from the Sun. The effective temperature and metallicity 
are T,ff = 5921 + 26 and [Fe/H] = 0.00 + 0.03. For this star 
we derived a mass M» = 1 .04 + 0.03 M© with a luminosity of 
1.24 Lq, an age of 3.5 + 1.9 Gyr, and a radius of 1.07±0.04. The 
log/?^^ of -4.98 ± 0.05 derived from CORALIE spectra leads 
to a calibrated stellar rotation period of Prot = 20 + 3 days. 



3. Radial velocities and orbital solutions 

The CORALIE echelle spectrograph was mounted on the 1.2 m 
Euler Swiss telescope in June 1998. It went through a major 
hardware upgrade in June 2007 to increase overall efficiency. 



The reader is referred to Segransan et al. 



( 2009 1 for more details. 



This hardware modification has affected the instrumental radial 
velocity zero point with offsets that depend on the target spec- 
tral type. For this reason, we decided to consider the upgraded 
spectrograph as a new instrument called CORALIE-07 (C07). 
We refer to the original CORALIE as CORALIE-98 (C98), and 
we adjust an offset between radial-velocities gathered before and 
after June 2007. Time series velocity data are fit with a Keplerian 
model using a genetic algorithm tool developed at the University 
of Geneva Observatory and called Yorbit ( Segransan et al. 201 1|. 
The interest of the genetic approach is to explore the different 
parameter spaces by creating a population of candidates to the 
Keplerian model and working toward better solutions. This tool 
has been developed and is especially useful for multiplanet sys- 
tems, but it delivers accurate results in single-planet systems. 
This technique avoids the fitting algorithm becoming trapped in 
a local minimum. The best individuals of the evolved popula- 
tion are then used as input parameter to a Levenberg-Marquardt 
algorithm to obtain the best minimization fit between the as- 
sumed Keplerian model and the observed data. The parameter 
confidence intervals computed for a 68.3% confidence level af- 
ter 10000 Monte Carlo trials define the uncertainties. A Markov 
chain Monte Carlo (MCMC) algorithm is used to estimate the 
uncertainties of the fitted parameters when the orbital period is 
longer than the time span of the observations. 
Because some of the planetary companions presented in this pa- 
per are fairly massive, the knowledge of the radial-velocity pa- 
rameters can be used to search for a Hipparcos astrometric sig- 
nature in the same way as in Sahlmann et al. (2011 1. Except for 
HD 10647b (see Sect. 4.1 1 Hipparcos astrometry was not able to 
constrain the upper mass of the planets. For the lightest plan- 
etary companions, the astrometric signal is probably too weak, 
while for the longest orbital period, the comparatively short time 
span of the Hipparcos data implies poor orbit coverage (less than 
25%). The standard orbital parameters derived from the best fit- 
ted Keplerian solution to the data are listed in Table|2] 

3.1. A six-year period Jovian planet around HD 27631 

We have acquired 60 radial velocities measurements of 
HD 27631 since November 1999. Twenty-three of these were 
obtained with CORALIE-98 with a typical signal-to-noise ra- 
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Table 2. Single-planet Keplerian orbital solutions for HD 2763 1 , HD 98649, HD 1 065 1 5 A, HD 1 66724, HD 1 96067, and HD 220689 
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tio of 40 (per pixel at 550 nm) leading to a mean measurement 
uncertainty including y-noise and calibration errors of 7.8 ms"'. 
Thirty-seven additional radial-velocity measurements were ob- 
tained with CORALIE-07 with a mean signal-to-noise ratio of 
75 leading to a mean measurement uncertainty of 6.1ms-'. 
Figure|3] shows the CORALIE radial velocities and the corre- 
sponding best-fit Kepleiian model. Since each data set has its 
own velocity zero point, the velocity offset between CORALIE- 
98 and CORALIE-07 is an additional free parameter to include 
in the model, as described before. The resulting orbital param- 
eters are P = 2208±66 days, e = 0.12±0.06, and K = 23.7±1.9 
m-', implying a minimum mass nip sin / - 1.45±0. 14 Mjup or- 
biting with a semi-major axis a - 3.25 ± 0.07 AU. The rms to 
the Keplerian fit is 6.31ms-' for CORALIE-98 and 7.56 ms-' 
for CORALIE-07, yielding a reduced of 1.37 ±0.23 and 1.77 
as goodness of fit. The single-planet Keplerian model thus ade- 
quately explains the radial-velocity variation, because for each 
instrument the rms is compatible with the mean measurement 
uncertainties. 



3.2. A 7 Jupiter-mass companion on an eccentric orbit of 14 
years around HD 98649 

Eleven radial-velocity measurements of HD 98649 were ob- 
tained between February 2003 and March 2007 with CORALIE- 
98 showing a slight linear drift of -9.7 ms-'yi"-'. The average 
measurement unceitainty is 7.7 ms ' for a mean signal-to-noise 
ratio of 48. Since the trend was very constant and the radial- 
velocity dispersion did not show any sign of a potential plan- 
etary companion, the observing frequency was decreased. In 
December 2008 after nearly six years of observation, a new mea- 
surement made with CORALIE-07 revealed a radial- velocity 70 
ms-' below the expected value. 

Additional radial-velocities gathered with CORALIE-07 un- 
veiled a companion with a minimum mass ;«pSin/ - 6.8±0.5 
Mjup on an eccentric orbit e - 0.85±0.05 around HD 98649 



with a period of P = 4951!;4gj days, implying a semi-major axis 
a = 5.6 + 0.4 AU. The thirty-two radial-velocities gathered with 
CORALIE-07 have a typical signal-to-noise ratio of 85, lead- 
ing to a mean measurement uncertainty of 5.9 ms-' . With an ob- 
serving time span of 9.4 years the radial velocity measurements 
cover only 70% of orbital phase, explaining why the period is 
poorly constrained. However, the good data coverage near the 
periastron passage of the planet allows us to constrain the orbital 
eccentricity and the velocity semi-amplitude very well. Figure 
|5]shows the covariance between eccentricity, semi-amplitude K, 
and orbital period for 10000 Monte Carlo realizations. 
An increased eccentricity can lead to a longest orbital period 
with a slightly higher semi-amplitude. However, the lowest limit 
of these two parameters is quite sharp and a more circular or- 
bit does not change K and P. With less than an orbital phase 
coverage, we can not exclude a long-term trend in the mea- 
surements. Adding a positive linear drift to the model leads to 
a slight decrease in the peiiod (down to a factor of 0.8) but 
does not change the eccentricity which stays above e < 0.82. 
Applying the Akaike information criterion ([Akaike 1974 1, which 
attempts to find the model that best explains the data with a 
minimum of free parameters, reveals that the single Keplerian 
model without drift fits our measurements better (Fig|4]l. The 
residuals around the single planet orbital solution for both data 



sets io-(o- 



C)c98 



10.31 ms , o-(o-c)c 



5.68 ms ) agree 



with the precision of the measurements. With a reduced of 
1.35±0.27 and a G.o.F = 1.42, the single-planet Keplerian model 
adequately explains the radial-velocity variation. HD 98649b is 
in the top five of the most eccentric planetary orbits and the most 
eccentric planet known with a period longer than 600 days. In 
the absence of detection of an outer companion, we cannot in- 
voke the Kozai mechanism (Koz ai]|T962| l as an explanation for 
the origin of the high eccentricity of the orbit. With an orbital 
separation of 10.4 AU at apoastron, leading to an orbital sep- 
aration up to 250 milliarcsecond, HD 98649b is an interesting 
candidate for direct imaging. 
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Fig. 3. Radial-velocity measurements as a function of Julian 
Date obtained with CORALIE-98 (red) and CORALIE-07 (blue) 
for HD 2763 1 . The best single-planet Keplerian model is repre- 
sented as a black curve. The residuals are displayed at the bottom 
of the top figure and the phase-folded radial-velocity measure- 
ments are displayed in the bottom diagram. 



Fig. 4. Radial-velocity measurements as a function of Julian 
Date obtained with CORALIE-98 (red) and CORALIE-07 (blue) 
for HD 98649. The best single-planet Keplerian model is repre- 
sented as a black curve. The residuals are displayed at the bottom 
of the top figure and the phase-folded radial-velocity measure- 
ments are displayed in the bottom diagram. 



3.3. A ten Jupiter-mass planet on a ten years orbit around 
HD 10651 5 A 

HD106515A has been observed with CORALIE at La Silla 
Observatory since May 1999. Twenty-two Doppler measure- 
ments were gathered with CORALIE-98 with a typical signal- 
to-noise ratio of 47 (per pixel at 550 nm), leading to a mean 
measurement uncertainty of 8.4 ms"'. Twenty-four additional 
radial-velocity measurements were obtained with CORALIE- 
07 with a mean signal-to-noise ratio of 61 leading to a mean 
measurement uncertainty of 6.3 ms ^ HD 1065 15Ab has com- 
pleted slightly more than one orbital revolution since our first 
measurement (Fig|6]l, and the phase coverage near the perias- 
tron has been greatly improved by the recently gathered obser- 
vations. Thus our Monte Carlo realizations for a single-planet 
Keplerian model point to well constrained orbital elements for 
HD 1065 15Ab. The Icr confidence interval for K is 158.2±2.6 
ms~', which leads to a planetary mass of 9.61 ±0.14 Mjup. For the 
same confidence interval, the eccentricity and period distribution 



results in the values e = 0.572+0.011 and P = 3630+12 days 
with a semi-major axis of 4.590+0.010 AU. The orbital separa- 
tion ranges from 1.96 AU at periastron to 7.22 AU at apoastron 
which corresponds respectively to 56 and 205 milliarcsecond. 
The residuals to the fit are somewhat greater than the precision 
of the measurements (o-^o-c)cm - 9-63 ms"', cr^o-Ocm - 6-66 
ms ') and are probably caused by a combination of stellar jitter 
and light contamination of the second component of the binary. 
However, with a reduced of 1.55+0.28 and a G.o.F = 2.16, the 
data are well modeled with a single Keplerian orbit. The pres- 
ence of a stellar companion suggests that the eccentric orbit of 
HD 106515Ab may arise from a Kozai mechanism (see Sect.|5]|. 

3.4. HD 166724b: a 1 4-year period companion of 3.5 Jupiter 
mass 

We began observing HD 166724 in September 2001 and have 
acquired 99 Doppler measurements for this target since then. 
Eighteen of these were obtained with CORALIE-98 with a 
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Fig. 5. Covariance between orbital period, semi-amplitude, and 
eccentricity for HD 98649b (top), HD 166724b (middle), and 
HD 219077b (down) as characterized after 5x10^ Markov-chain 
Monte Carlo iterations 



signal-to-noise of 27 leading to a mean measurement uncer- 
tainty of 10.6 ms"'. Twenty-nine more radial velocities were 
gathered with CORALIE-07 and have a mean final error of 7.3 
ms"' with S/N = 51 per pixel at 550 nm. Since HD 166724 is 
also part of the HARPS high-precision survey searching for low- 
mass planets around dwarf stars, fifty-five radial-velocities with 
higher signal-to-noise and precision have been obtained with this 
spectrograph since April 2004. The HARPS data on HD 166724 
have a mean measurement uncertainty (including photon noise 
and calibration errors) of 0.95 ms"' corresponding to a typical 
signal-to-noise ratio per pixel at 550 nm of 91. To account for 
the slightly high stellar activity (log/?^^ = -4.73), a velocity 
jitter of 3.5 ms ' was quadratically added to the mean radial ve- 
locity uncertainty. 

Even though the observing time span of 1 1 years does not cover 
the putative complete orbital period, the shape of the orbit is well 
constrained by the Doppler measurements gathered near and af- 
ter the passage through periastron, which constrain the eccentric- 
ity and the semi-amplitude. However, a second observation at the 
periastron passage is needed to definitively constrain the orbital 
period and the other orbital parameters. Running Yorbit followed 
by MCMC simulations led to a period of 14. l;^[ j years with an 
eccentricity of e = 0.734±0.020. A correlation (R=0.7) remains 
between eccentricity and period (Fig|5]), however the value of 
semi-amplitude {K - 71.0+1.7 m"') is decorrelated from these 
two parameters. The corresponding best Keplerian fit and time 
series radial velocities using C98, C07, and HARPS are plotted 
in FigjT] The single-planet Keplerian fit yields a reduced of 
1 .24±0. 16, however the residuals around the model remain fairly 
high with an rms = 3.57 ms"' for HARPS. The large re-emission 
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Fig. 6. Radial-velocity measurements as a function of Julian 
Date obtained with CORALIE-98 (red) and CORALIE-07 (blue) 
for HD 1065 15A. The best single-planet Keplerian model is rep- 
resented as a black curve. The residuals are displayed at the bot- 
tom of the top figure and the phase-folded radial-velocity mea- 
surements are displayed in the bottom diagram. 



peak in the Call region and the relatively high value of log R'^^ 
observed for HD 166724 point toward influence of the stellar jit- 
ter to explain this remaining excess of variability. 

3.5. HD 196067b: a Jovian companion (6.9 Mj^p) on a 
ten-year orbit 

Eighty-two observations have been obtained on HD 196067 
since September 1999. Thirty radial velocities were gathered 
with CORALIE-98 and fifty-two with CORALIE-07. The mean 
signal-to-noise ratio and measurement uncertainty are 72 and 
7.5ms-i for C98 and 115 and 5.8 ms-' for C07. The time se- 
ries velocities and the corresponding best-fit Keplerian orbit are 
plotted in Figjs] The rms to the Keplerian fit is 8.88 ms-' for 
C98 and 9.92 ms ' for C07 with G.o.F = 6.55 and a reduced 
of 2.48±0.26. The residuals around the single Keplerian model 
are slightly more than the precision of the measurements for 
C07, which points toward the stellar jitter or light contamina- 
tion and/or guiding errors induced by the second component of 
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Fig. 7. Radial-velocity measurements as a function of Julian 
Date obtained with CORALIE-98 (red), CORALIE-07 (blue), 
and HARPS (yellow) for HD 166724. The best single-planet 
Keplerian model is represented as a black curve. The residuals 
are displayed at the bottom of the top figure and the phase-folded 
radial-velocity measurements are displayed in the bottom dia- 
gram. 



Fig. 8. Radial-velocity measurements as a function of Julian 
Date obtained with CORALIE-98 (red) and CORALIE-07 (blue) 
for HD 196067. The best single-planet Keplerian model is repre- 
sented as a black curve. The residuals are displayed at the bottom 
of the top figure and the phase-folded radial-velocity measure- 
ments are displayed in the bottom diagram. 



the binary. Due to the poor coverage during the closest approach, 
the data are somewhat overfitted for C98. This explains why the 
rms to the fit is better for this instrument. 

The radial velocity of HD 196067 remained very stable during 
the first three years of measurement and the observing frequency 
of the target was diminished. This explains the poor data sam- 
pling before the periastron. Because HD 196067b has completed 
one orbital revolution since the beginning of the CORALIE sur- 
vey, we have only gathered two observations near the radial ve- 
locity minimum. This therefore allows us to compute the semi- 
amplitude with a 68.3% confidence interval of 80 to 256 ms"', 
yielding to a derived planetary mass of 5.8 to 10.8 Mjup. The 
period is, however, not perfectly constrained with a minimum 
value of 9.5 years corresponding to an eccentricity of e=0.57. 
The choice of a higher eccentricity of 0.84 tends to increase the 
period up to 10.6 years. Finally, the Kozai mechanism induced 
by the second stellar component of the system can be invoked to 
explain the eccentricity of the planetary orbit (see Sect.|5]l. 



3.6. HD 219077b: an ten Jupiter-mass planet on a 
fifteen-year orbit 

HD 219077 has been observed with CORALIE at La Silla 
Observatory since June 1999. Ninety-three Doppler measure- 
ments were gathered on this target during the past 13.3 years 
with C98, C07, and HARPS. The typical signal-to-noise ratio 
(per pixel at 550 nm) for each instrument is respectively 46, 112, 
and 211 corresponding to a mean measurement uncertainty of 
8.3,5.7,1.0 ms"'. As HARPS measurements are contemporary 
with both CORALIE data sets, the velocity offsets between the 
spectrographs are therefore well constrained. The radial veloc- 
ities time series display a linear drift of 13 ms^'yr"' over the 
first eight years of observation. This trend was compatible with 
a brown dwarf companion at a distance that would make the de- 
tection possible with the resolution of the Very Large Telescope. 
A direct observation of the companion of HD 219077 was at- 
tempted in 2006 by Montagnier et al. (in prep) with the simul- 
taneous differential imaging mode of the adaptive optics instru- 
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Fig. 9. Radial-velocity measurements as a function of Julian 
Date obtained with CORALIE-98 (red), CORALIE-07 (blue) 
and HARPS (yellow) for HD 219077. The best single-planet 
Keplerian model is represented as a black curve. The residuals 
are displayed at the bottom of the top figure and the phase-folded 
radial-velocity measurements are displayed in the bottom dia- 
gram. 



Fig. 10. Radial-velocity measurements as a function of Julian 
Date obtained with CORALIE-98 (red) and CORALIE-07 (blue) 
for HD 220689. The best single-planet Keplerian model is repre- 
sented as a black curve. The residuals are displayed at the bottom 
of the top figure and the phase-folded radial-velocity measure- 
ments are displayed in the bottom diagram. 



ment NACO. No detection was reported in their high-contrast 
images, and they can exclude a companion that is more massive 
than 100 Mjup. Combining their observations and the linear drift 
in the CORALIE data, they predict that HD 219077b has a prob- 
ability of 0.75 to be a brown dwarf. 

Since 2010, an acceleration in the radial velocity was observed in 
the time series of CORALIE-07. In September 201 1, the cadence 
of Doppler measurements was increased in order to optimize the 
data sampling near the periastron passage. Even if the orbital 
period is still not covered today, this strategy has allowed the or- 
bital parameters to be constrained through MCMC simulations. 
The C98, C07, and HARPS data sets and the corresponding best 
Keplerian fit are plotted in Fig.|9] The remaining correlations be- 
tween eccentricity, period, and semi-amplitude, characterized by 
the MCMC method, are presented in Fig.|5] Even if the semi- 
amplitude is well constrained (K = 181.4+0.8 m"') and the 1-cr 
confidence interval of the eccentricity is relatively small (0.767 
up to 0.773), the corresponding period range remains around 



0.7 years (14.7 yr < P < 15.4 yr), implying a minimum mass 
nip sin i for HD 219077b between 10.30 and 10.48 Mjup. 

3. 7. HD 220689b: a Jupiter analogue on a six years orbit 

The first CORALIE measurement of HD 220689 was gathered 
in September 2000. The first eleven of these were obtained with 
CORALIE-98 with a signal-to-noise of 52 and a mean measure- 
ment uncertainty of 8. 1 ms"' . Thirty-three more radial velocities 
were gathered with CORALIE-07 and have a mean final error of 
5.8 ms"^ with S/N = 99 per pixel at 550 nm. A very significant 
peak with FAP < 10"^ is identified around 2200 days in the 
Lomb Scargle periodograme ( GiUiland & Baliunas|1 987) of the 
CORALIE time series. This signal corresponds to a Jupiter-mass 
planet (m^ sin/= 1.06±0.09 Mjup) orbiting HD 220689 with P = 
2209+^f days on a sHghtly eccentric orbit (e = 0.16+[},'°). The 
residuals to the Keplerian fit are 6.08 ms ' for the combined data 
set leading to a reduced of 1.14+0.25. The semi-amplitude K 
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Table 3. Observed and inferred stellar parameters for HD 10647, 
HD 30562, and HD 86226. 



Parameters 




HD 10647 


HD30562 


HD86226 


Spectral type*"' 




F9 V 


F8 V 


G2Y 


lAci) 




5.52 


5.77 


7.93 


B — 




0.551 


0.631 


0.647 


7T' ' 


[mas] 


57.36+0.25 


37.84+0.35 


22.2+0.8 


My 




/I Q 1 

4. ji 


3.O0 


/I AA 


^eff 














4.62±0.04<''' 


4.34+0. 10<^' 


4.36+0.04*''' 


[Fe/H] 


[dex] 


O.OO+O.Ol'* 


0.29+0.06*''' 


0.05+0.03*''' 


vsin (0*"* 


[kms-'] 


4.9 


4.1 


24 




[Mo] 


1.11 ±0.02 


1.23+0.03 


1.06+0.03 


Lis) 


[LoJ 


141 


2.63 


1.06 


R(f) 


[Rol 


1.10±0.02 


1.53+0.04 


1.02+0.03 






-4.77±0.01<''> 


-5.05+0.09''' 


-4.90+0.05''' 


p(i) 


[days] 


10±3 


26+4 


23+4 




[Gyr] 


1.4±0.9 


3.6+0.5 


1.5+1.3 



Notes. <"' Parameter from HIPPARCOS (ESA'l997V <*' Parallax from 
the new Hipparcos reduction (van Leeuwen 20071. ''' Parameter de- 
rived us ing CORALIE spec tra. *"' Parameter derived using HARPS 
spectra jSousa et al.||2008f . *"' Parameter derived using CORALIE 
CCF (Santos et al. 20011. Parameter derived from Girardi et al. 
2000) models. *-''' The bolometric correction is com puted from Flower 
1996^ . **' Parameter derived using HARPS spectra ( |Lovis et al. 201 1 1. 
From the calibration of the rotational period vs. activity ( |Mamajek &| 
|Hillenbrand|2008l l. 



= 16.4+1.5 m~' of HD 220689b is one of the weakest radial ve- 
locities signals unveiled with CORALIE until today. However, 
the detection of the companion remains very robust with a sig- 
nal semi-amplitude around three times bigger than the rms to the 
model. No more significant power is detected in the time series 
periodograme after subtracting HD220689b's orbit. 

4. Updated parameters for three known exoplanets 

4.1. HD 10647Ab (HIP 7978b), a Jupiter-mass planet on a 
thousand-day period 

HD 10647b was first announced in 2003 by Mayor et al. in the 
XlVth lAP Colloquium. At this time, the small semi-amplitude 
of the orbit and the relatively high residuals to the Keplerian 
fit due to stellar activity, pushed the authors to carefully gather 
more data before publication. Three years later, ' Butler et al.| 
( 2006| l published orbital parameters for HD 10647b based on 
28 radial velocities gathered with the UCLES on the Anglo- 
Australian Telescope. However, the CORALIE radial veloci- 
ties and parameters were never published before this paper 
HD 10647 is a bright F9 dwarf located in the Eridanus con- 
stellation at 17.43 + 0.08 pc from the Sun. The star properties 
are summaiized in Table[3] Even though no clear prominent re- 
emission feature in Ca II absorption region can be seen in the 
HD 10647 spectra (Fig.|2]l, the logR'^j^, value of -4.77 ± 0.01 
points to moderate chromospherical activity. The possibly in- 
duced stellar spots combined with a slightly high measured pro- 
jected rotational velocity of vsin(i) - 4.9 kms"' led to ex- 
pecting some radial-velocity jitter on a typical timescale on 
the order of the rotational period = 10 + 3 days). A 

clear signature around P=1000days (Fig{TT]a) is identified in 
the Lomb Scargle periodogram of the 108 Doppler measure- 
ments gathered with CORALIE since September 1999. This sig- 
nal corresponds to a Jupiter-mass planet (irip sin / = 0.94±0.08 
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Fig. 11. Lomb-Scargle periodogram of the CORALIE radial ve- 
locities for HD 10647 (a). The dotted line indicates the 10"^ false 
alarm probability {3cr) and the dotted line the coiTesponding Icr 
limit. The peak at 267 days, is the one-year alias of the plane- 
tary signal at P=1000 days. It disappears in the residual's Lomb- 
Scargle (b), which no longer shows any signal. 



Mjup) on a 989.2±8.1 days orbit with a semi-major axis of 



=2.015+0.01 1 AU. Figure[T2|shows the combined C98 and C07 
times series with the corresponding best-fit Kepleiian model and 
the orbital radial velocity signature, as expected from the pa- 
rameters published by Butler et aL] ( |2006 1. The resulting orbital 
elements for HD 10647b are listed in Table|4] The residuals to 
the adjusted single-planet Keplerian model show a level of vari- 
ation that is slightly larger than expected from the precision of 
the measurements. This points toward the influence of the stellar 
jitter since the periodogram of the residuals (Fig{TT]b) no longer 
shows any signal with significant power Additionally, a con- 
straint on the radial-velocity solution can be provided by the 
intermediate astrometric data of the new Hipparcos reduction 
( van Leeuwen||2007^ . With the 123 measurements gathered by 
Hipparcos over 1128 days (1.14 orbits), we applied the method 
described by 'Sahl mann et all] ( |2011) l to search for the astromet- 
ric signature. Hipparcos can constrain the maximum mass of the 
companion of HD 10647 to 12.33 Mjup and thus exclude a brown 
dwarf or M dwarf companion. 



4.2. Independent confirmation and refined orbital parameters 
of the long-period and eccentric Jovian planet 
HD 30562b 

HD 30562 is identified as a bright F8V in the Hipparcos cat- 
alog with an apparent magnitude of V = 5.77 at a distance 
of 26.43 ± 0.24 pc from the Sun. The spectroscopic analysis 
of the CORALIE spectra results in an effective temperature 
T,ff = 5994 + 46 and metallicity of [Fe/H] = 0.29 + 0.06. For 
this star we derived a mass M» - 1 .23 + 0.03 Mg with an age of 



3.6 ± 0.5 Gyr. |HaII et ar| ( [2007l ) find that HD 30562 is chromo- 
spherically inactive with an activity index of logR'^^, - -5.15 
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Fig. 12. Radial-velocity measurements as a function of Julian 
Date obtained with CORALIE-98 (red) and CORALIE-07 (blue) 
for HD 10647. The best single-planet Keplerian model is repre- 
sented as a black curve while the red dotted line represents the 
orbital radial velocity signature as expected from the parameters 
published by Butler et al.| (|2006). The residuals are displayed at 
the bottom of the top figure and the phase-folded radial- velocity 
measurements are displayed in the bottom diagram. 



Fig. 13. Radial-velocity measurements as a function of Julian 
Date obtained with CORALIE-98 (red) and CORALIE-07 (blue) 
for HD 30562. Data from Lick Observatory Fischer et al. (2009| 
are plotted in yellow. The best single-planet Keplerian model 
is represented as a black curve. The residuals are displayed at 
the bottom of the top figure and the phase-folded radial-velocity 
measurements are displayed in the bottom diagram. 



in agreement with the value of log 7?^^ - -5.05 + 0.09 derived 
in the present paper from CORALIE spectra. We began observ- 
ing HD 30562 in January 1999 and have acquired 130 Doppler 
measurements for this target since then. The first twenty-two of 
these were obtained with CORALIE-98 with a signal-to-noise 
of 68 leading and a mean measurement uncertainty of 7.1 ms"'. 
Sixty -four more radial velocities were gathered with CORALIE- 
07 and have a mean final error of 5.8 ms"' with S/N =130 per 
pixel at 550 nm. In August 2009, HD 30562b was announced 
by Fischer et al. (2009) based on forty-four data points gath- 
ered at the Lick Observatory. Since a good phase coverage is 
needed to constrained the high eccentricity of the orbit, we de- 
cided to combine the three data sets to refine the orbital param- 
eters. The corresponding best Keplerian fit and time series ra- 



dial velocities from C98, C07, and Lick are plotted in Fig. 13 
Since each instrument has its own velocity zero point, the ve- 
locity offsets between the three spectrographs are additional 



free parameters in the model. The data sets have contempo- 
rary observations, so the velocity offsets are well constrained. 
The residuals to the adjusted single-planet Keplerian model 
show a level of variation of 0--6J6 ms"', yielding a reduced 
of 1.71±0. 17. The rms to the Keplerian fit is 4.97 ms"' for 
CORALIE-98, 6.88 ms ' for CORALIE-07, and 7.42 ms-' for 
Lick data. The resulting orbital parameters are f = 1 159.2±2.8 
days, e = 0.778+0.013, K = 36.8±1.1 m"', implying a mini- 
mum mass nip sin i - 1.373+0.047 Mjup orbiting with a semi- 
major axis 0=2.315+0.004 AU. These results are consistent with 
the parameters by [Fischer et al. ( |2009 1. The combined fitting 
of the data sets leads to a sUghtly higher eccentricity and semi- 
amplitude, implying a planetary mass around 5% above the value 
previously announced. Interestingly, due to the high eccentricity, 
the transit probability of HD 30562b is quite high (x!l%) for a 
planet with a revolution period in this range. 
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Table 4. Single planet Keplerian orbital solutions for HD 10647, 
HD 30562, and HD 86226. 
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Fig. 14. Radial-velocity measurements as a function of Julian 
Date obtained with CORALIE-98 (red) and CORALIE-07 (blue) 
for HD 86226. The best single-planet Keplerian model is repre- 
sented as a black curve while the red dotted line represents the 
orbital radial velocity signature as expected from the parameters 
published by Min niti et al.|p009 1. The residuals are displayed at 
the bottom of the top figure and the phase-folded radial- velocity 
measurements are displayed in the bottom diagram. 



Parameters 




HD 10647 


HD30562 


HD86226 


F 


[days] 


989.2±8.1 


1159.2±2.8 


1695±58 


K 


[m-'] 


18.1 + 1.7 


36.8±1.1 


15.3±1.7 


e 




0.15±0.08 


0.778±0.013 


0.15±0.09 


w 


[deg] 


212±39 


78.4±3.0 


193±61 


To 


[JD] 


2453654±99 


2453601. 4±2.8 


2454821 ±288 


a I sin i 


[10--' AU] 


1.71±0.15 


2.47±0.08 


2.49±0.28 


fiiin) 


[10-9 


0.70±0.18 


1.49±0.15 


0.73±0.24 


nip sin i 


[Mj„p] 


0.94±0.08 


1.373±0.047 


0.92±0.10 


a 


[AU] 


2.015±0.011 


2.315±0.004 


2.84±0.06 


Jem 


[kms-'] 


27.6782 


77.2364 


19.7389 






±0.0011 


±0.0021 


±0.0039 


AVc07-C98 


[ms-'] 


-2.4±2.2 


2.9±24 


-3.0±4.0 



AVifct, 

MIKE 



[ms->] 



-77237. 8±2.3 -19734.4±4.1 



N 




108 (92/16) 


130 (22/64/44) 


65 (11/41/13) 


AT 


[years] 


13.0 


13.7 


13.0 






1.65±0.18 


1.71±0.17 


1.43±0.22 


G.o.F 




3.90 


4.64 


2.09 


O" (O-C) 


[ms-'] 


8.92 


6.76 


6.88 




[ms-i] 


9.38 


4.97 


8.70 


cr (O-C)coi 


[ms-i] 


7.51 


6.88 


6.88 


^'MIKE 


[ms-'] 




7.42 


5.01 



Notes. Confidence intervals are computed for a 68.3% confidence level 
after 10000 Monte Carlo trials. C98 stands for CORALIE-98 and C07 
for CORALIE-07. AT is the time interval between the first and last 
measurements, is the reduced G.o.F. is the Goodness of Fit and 
<T(o-C) the weighted rms of the residuals around the derived solution. 



dial velocity decrease predicted by published orbital parameters. 



Figure 14 shows the combined C98, C07, and MIKE times series 
with the corresponding best-fit Keplerian model and the orbital 
radial velocity signature, as expected from the parameters pub- 
lished by Minniti et al.| ( i2009 ). The new CORALIE observations 
confirm the presence of a Jovian planet {nip sin i - 0.92+0.10) 
on a 4.6-year period around HD 86226. However, the combined 
fitting of the data sets leads to an orbit eccentricity and a much 
lower semi-amplitude than previously announced. The resulting 
orbital elements for HD 86226b are listed in Table|4] The case 
of HD 86226b is a good illustration of the bias toward overes- 
timated Mp sin / and eccentricity for a moderate number of ob- 
servations and relatively low signal-to-noise ratio of K/cr < 3 
dShen & Tumer|2008i l. 



4.3. Independent confirmation and new orbital parameters of 
the long-period Jupiter-mass planet around HD 86226 

HD 86226 has a spectral type G2V with an astrometric paral- 
lax of ;r = 22.2 + 0.8 mas, which sets the star at a distance 
of 45.0 + L6 pc from the Sun in the Hydra constellation. The 
star is chromospherically quiet with a log 7?^^ of -4.90 ± 0.05 
leading to a stellar rotation period of Prot = 23 ± 4 days. The 
first measurement of HD 86226 was gathered in March 1999. 
Thirteen observations were obtained with CORALIE before the 
announcement of HD 86226b made by Min niti et al.| ( |2009| l 
based on thirteen MIKE measurements. At this time, despite a 
small excess of variability, no significant signal could be un- 
veiled in the CORALIE data. Since March 2010 a new ob- 
servation campaign has been carried out on HD 86226, lead- 
ing to thirty-nine new Doppler measurements. The observing 
frequency was relatively high in order to observe the fast ra- 



5. Discussion & conclusions 

We have reported in this paper the discovery of seven giant plan- 
ets with long orbiting periods discovered with the CORALIE 
echelle spectrograph mounted on the 1.2-m Euler Swiss tele- 
scope located at La Silla Observatory. Interestingly, we notice 
that the hosts of these long-period planets show no metallicity 
excess. Five of them - HD 98649b, HD 106515Ab, HD 166724b, 
HD 196067b, and HD 2 19077b - have periods around or above 
ten years, and they constitute a substantial addition to the pre- 
viously known number of planets with periods in this range. 
These five exoplanets have surprisingly high eccentricities above 
e > 0.57 (Fig.[T5]). HD 98649b, HD 166724b, and HD 219077b 
(e - 0.85, e = 0.734, and e - 0.770) are even the three most ec- 
centric planets with a period longer than five years. Our Monte 
Carlo and MCMC simulations show that this parameter is well 
constrained for these targets and that no bias toward overesti- 
mated eccentricities ( |Shen & Tumer|2008] l could be invoked. We 
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Fig. 15. Observed orbital eccentricity vs. period for the known 
extra-solar planets from the exoplanet.eu database (Schneider 
et al.pOlT I. The parameters of the seven companions unveiled 
in the present paper are plotted in red. 



notice that, due to their high eccentricity orbits, the orbital sepa- 
ration range up to 10.4 and 1 1 .0 AU at apoastron for HD 98649b 
and HD 219077b. This can lead to an angular separation up to 
250 and 375 milliarcsecond, which makes two interesting can- 
didates for direct imaging. In both host stars HD 1065 15A and 
HD 196067, a third massive body is present in the system. To in- 
vestigate whether the eccentricity of the orbit may arise from the 
Kozai mechanism ( Kozai|1962|), we estimate t heir Kozai period 
of eccentricity oscillation ( |Mazeh & Shaham||1979| |Ford et aT] 
[2000) 1 



Kozai 



plane! 



(1) 



where a and PpUmei are the semi-major axis and the period of 
the planetary orbit, Mp and M, the mass of the primary and the 
secondary stars, ai, and the semi-major axis and the eccentric- 
ity of the binary. Assuming as eccentricity distribution /(e) - 2e 
for the binaries with P > 1000 days ( [Duquennoy & Mayor] 1 99 ij l 
and taking the median of the distribution (e^, - ^), we obtain 

PKozai = 3.3 X 10'' years for HD 1065 15Ab and PkozoI = 3.5 x 10^ 
years for HD 196067. The period of eccentricity oscillations 
PKozai has to be compared to age of the system and to the apsi- 
dal precession due to the relativistic correction to the Newtonian 
equation of motion ( |Holman et al.|199"7] l, which is given by 



rei 



planet ' 



3GM,,„r 



(2) 



This yields P^ei = 1.1 x 10^ years for HD 106515Ab and P„; = 
7.4 X 10** years for HD 196067b. For the two planets the eccen- 
tricity oscillation half-period is shorter than the age of the sys- 
tem and at least one order faster than the precession due to the 
relativistic effects. Since both binary stars have semi-major axis 



wider than 40 AU, we can assume that the planet and the bi- 
nary orbital planes have uncorrelated inclinations ('Halel|1994 1. 
Thus the Kozai pumping mechanism is a viable explanation for 
the eccentricities of the HD 106515Ab and HD 196067b orbits. 
Because the presence of a third massive body cannot be inferred 
from the data of HD 98649b, HD 166724b, and HD 219077b, the 
origin of the eccentricity of these systems remains unknown. 

Interestingly, half of the known planets with a semi-major axis 
larger than 4 AU (15/30) are part of multiple systems, whereas 
this ratio is much lower (~25%) for the companions in the 1-4 
AU range. The most likely explanation for this is a detection bias 
induced by the follow-up of the first detection of a smaller pe- 
riod planet, which triggers a follow-up of the star and tends to 
increase the span of the time series. For many configurations of 
multiple systems, high eccentricities are incompatible with sta- 
bility, so these systems are uncommon. As a consequence, the 
small number of single companions beyond 4 AU biased the ec- 



centricity distribution toward low eccentricities (Fig. 16 1 



This trend is probably reinforced when one considers that very 
eccentric planets are easily overlooked. If one observes the flat 
part of an eccentric orbit one may detect only small variations, 
and many surveys would cast the star aside as uninteresting. This 
point and the large amount of data gathered by the CORALIE 
survey do not explain the high fraction of very eccentric planets 
unveiled in the present paper. However, these are small number 
statistics, and a Kolmogorov-Smirnov test (Massey 1951| ) ap- 
plied to the eccentricity distribution in the range 1-4 AU, and 
beyond, gives a good probability (~50%) for the null hypothesis 
(i.e. that the samples are drawn from the same distribution) with 
or without the planets presented in this paper. 

The same remark applies to the lack of metallicity excess of the 
hosting stars (of this paper). Even though Boisse et al. (2012 1 
has noticed that beyond 4 AU, the well described greater occur- 
rence of gas giants around metal-rich stars ( [Santos et al.||2004| 
[Fischer & Valenti 2005 ; Mayor et al. 201 1 1 seems to remain true. 



the current statistic is still very poor Interestingly, we notice in 
Fig.[T6)that this correlation appears to be clearer (below and be- 
yond 4 AU) for multiple systems than for single ones. However, 
an unbiased sample and a statistical approach taking account of 
the detection efficiency for each star would be needed to confirm 
this hypothesis. 

In addition to the seven new companions presented in this paper, 
we also have published refined orbits for three already known ex- 
oplanets - H D 10647b (Butler et al. 2 006|), H D 30562b ( Fischer[ 
|etaL][2009l l, and HD 86226b (Minniti et al.[[2009.) . CORALIE 
time series combined with the published data sets, indepen- 
dently confirms the existence of the three companions. In the 
case of HD 86226b, our derived orbital parameters are, however, 
far from the published ones with an eccentricity decreasing from 
e - 0.73 to e = 0.15 and period more than hundred days longer. 
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Fig. 16. Distribution of the host star metallicities (left) and or- 
bital eccentricities (right) for the gas giant planets (>100 Mjg). 
The upper histograms represent the planets with semi-major axis 
between 1 and 4 AU and the lower quadrant is for planets beyond 
4AU. The hatched areas plot the contributions of the multiple 
systems and the red color is used for the 5 companions (with 
a > 4 AU) unveiled in the present paper The gray arrows mark 
in every quadrant the median of the distribution for the multiple 
systems and the black ones are for single-planet systems. 
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